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Введение. Если круглая пластина нагружена по боковой поверхности, возможны потеря ее устойчивости и 
выпучивание. Данную проблему можно решить с учетом бифуркации. В этом случае пластина рассматривается 
как нелинейно-упругое тело. В задачах устойчивости нелинейной упругости важен выбор зависимости между 
напряжениями и деформациями. В ранних работах, посвященных этой проблеме, рассматривались простые 
законы состояния (конститутивные уравнения). В качестве примера можно привести материал «гармонического 
типа» Сенсенига. 

Материалы и методы. Для круглой пластины из материалов Мурнагана и Блейтца и Ко получены уравнения 
нейтрального равновесия. В предположении однородной начальной деформации рассмотрена задача 
устойчивости этой пластины. Строгие трехмерные уравнения нейтрального равновесия позволяют исследовать 
смежные формы равновесия с учетом физической и геометрической нелинейности. Вывод этих уравнений 
основан на теории наложения малой деформации на конечную. 

Результаты исследования. Продвижение в решении соответствующего векового уравнения (с нелинейным 
вхождением параметра) для практически важных законов упругости Мурнагана и Блейтца и Ко возможно лишь 
с использованием численных методов. Разработанный метод расчета бифуркационных значений параметров 
нагрузки дает возможность проанализировать влияние нелинейности. 

Обсуждение и заключение. Исследовано влияние физической и геометрической нелинейности на величину 
верхнего критического значения параметра начальной деформации. Полученные результаты могут быть 
использованы при оценке достоверности модулей упругости третьего порядка для разных физических 
материалов. Данные об этих модулях пока малочисленны. Численные исследования показали, что к 
приведенным в некоторых источниках константам следует относиться с осторожностью. Обсуждается также 
использование модулей упругости в законе состояния Блейтца и Ко. 
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Введение. В настоящее время изучение новых сравнительно простых и адекватных законов состояния 
различных материалов, требующих учета нелинейных деформаций, представляет собой активно развивающееся 
направление механики сплошной среды [1-10]. В рамках теории наложения малой деформации на конечную 
выведены трехмерные уравнения нейтрального равновесия для материалов Мурнагана и Блейтца и Ко. На 
основе этих уравнений рассмотрен пример торцевого выпучивания равномерно сжатой круглой 
пластины [11, 12]. 

Материалы и методы 

Уравнения нейтрального равновесия. Пусть г, 9, 2 — цилиндрические координаты точки в 
недеформированном состоянии. Предполагаем, что начальная деформация тела определяется радиус-вектором 
В: 

В =А(ке, + 421, (1) 


где К(г) = аг ‚ параметры а, 4 — константы, е,, е» в — базисные векторы. 


Для координат этой точки в начально-деформированном состоянии имеем: 
В=аг, Ф=о Г=ф. 
Следовательно, градиент деформации, мера деформации Фингера и ее главные инварианты определены 
соотношениями: 
рт — ., — хр2 
УК=УЕ =а(е,е, +ерер) +4, Е =УК`, 


п =2а’+4*, Г =а+2а’4`, [=а\а?. (2) 
Здесь и далее У, У — набла-оператор в метрике недеформированного и начально-деформированного 


состояния: У =УЮ`".У. 


Будем использовать уравнения нейтрального равновесия, предложенные А. И. Лурье [13]: У@ =0, где 
тензор @ — линейный дифференциальный оператор над вектором добавочных перемещений У’. Выражения 
компонент этого тензора, конкретизированные с учетом законов состояния Мурнагана и Блейтца и Ко, 
получены в работах [14, 15]. В представлениях тензора © компоненты — некоторые функции, определяемые в 
результате решения краевой задачи начальной деформации. 

Результаты исследования. Для рассмотрения изгибной формы бифуркации равновесия пластины 
аналогично работам [11, 12] принимаем следующий вектор добавочного перемещения 

УИ! = и(г,2)е, + и(г,2)4 . (3) 
Учитывая (1), (2) и (3), тензор Лурье имеет вид: 
9 = Де, е, + Ве е, +Сы, + Се, + Нед. 
Здесь 
да а И вы Е Е Пао 
г д» 02 г д» 02 г д 





о 
02 д’ 02 Ге) 


Отметим, что структура этих операторов типична в задачах устойчивости цилиндрических нелинейно- 
упругих тел. Опуская выражения остальных коэффициентов, приведем, например, формулы для А. В случае 





закона состояния Мурнагана [13, 14] А выражается соотношением 


“(2 Зу, +4, 





А = 


У 
Е : четы 
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Здесь и далее ^. и и — модули упругости Ляме, у, у,, у, — константы упругости третьего порядка. 
Для материала Блейтца и Ко [15]: 
А и 2и(1 =. В) ра + В аа , 
ат В 


где В — уточняющий модуль упругости. 
Уравнение нейтрального равновесия равносильно системе дифференциальных уравнений относительно 


компонент вектора У’: 
2 1 2 2 
а [2 м [до вобу-о 








РИ а дубё а ' 02? 4 
2: 2 
ыы), о и 
02\ 0" г д гб’ а ` 02 


Предполагая, что по боковой поверхности пластина нагружена равномерным давлением, 
дифференциальные уравнения (4) дополним условиями равновесия на границе: 





и(",2)]_ =0, ©9,.|_, =0. 
На торцах, т. е. при 
р 
а 9.=0, 9. =0. (5) 
Применим подстановку 
и=Х, (2). Л (КГ), 6) 
и= 2, (2).Ло (А„г). 
Здесь п=1,2,..., /ь, /, — функции Бесселя нулевого и первого порядка, а числа К,’, — нули функции Л,(х). 
Примем х = Аг и применим равенства 
[9 4?.Л (К АЛ, (Е 
м ло, АО -юло, О ло --ьло, 
а’ Ч’ Ч’ 


а также известные тождества для функций Бесселя 
1 1 1 
п ' п" ' 
лы, ИЕЛыЕЛ, 
х х х 


Выполнив перечисленные выше действия, получим, что переменные в дифференциальных уравнениях 


(4) и граничных условиях (5) разделяются. Приходим к следующей краевой задаче для обыкновенных 
дифференциальных уравнений: 





А 
Хх, г г р К, С —_ 7, = 0, 
Н, С 
(7) 
[е; 6,+С, 


ЕЯ, +, 


п 





Х! =0. 
3 3 
СХ, (+#/2)+С,2, (+1/2)=0, (8) 
К,6,2, (+1/2)-С,Х! (+1/2) =0. 
Два однородных линейных обыкновенных дифференциальных уравнения (7) и четыре краевых 
условия (8) приводят к задаче на собственные значения с нелинейным вхождением параметра. В данной задаче 
таким параметром является величина (1 _ а) . Систему (7) приведем к нормальному виду: 


У =У,», 


А О. +С 
= уча, — у,, 
Н, \@р 





1 
Уз = У, 
С С, +С 
' 2 1 2 2 
= К, Уз К, У. 
С, С, 
Т Т 
Здесь приняты обозначения: аргумент 1=ЕФ, (х; у; У.) = (х,; ЖИ 2) . Пусть 
> Т 
фундаментальная система решений — это следующие четыре вектора у, =(у,;; у: У; У,). Тогда общее 
Т 4 
решение системы имеет вид (Х,; ЛХ’; 7,; 2) и = >, у, где & — произвольные постоянные. Зададим, 
= 


например, при 2=-й/2 начальные данные столбцами единичной матрицы четвертого порядка. Затем решим 
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численно задачу Коши с этими начальными условиями. В результате в правой части получим значения 
базисных функций, то есть векторы У;(+ й/2). Используя граничные условия (8), приходим к однородной 


системе линейных алгебраических уравнений относительно &, &,, 6, &: 
4 
а, =0, где 1=1, 2, 3, 4. 
= 


Здесь коэффициенты — элементы матрицы А: 
а, =ЮСу, (-1/2) + СУ. (-1/2) ‚а, = СУ (1/2) + Су (1/2) › 


а; =К, С, Уз; (-#/2)- СУ, (-1/2) ‚ 44, =К, СУ, (й/2) — б1у,; (1/2) у 
Однородная система имеет нетривиальное решение, если выполнено условие 
4е/А=0. (9) 
Выражение определителя включает параметры нагрузки а, 4, нули функции Бесселя, а также модули 


упругости Л, и, у, у,, ъ, (для материала Мурнагана) либо /., и, В (для материала Блейтца и Ко). 


12 


Параметры а, @, задающие начальную деформацию, взаимосвязаны. Осевая сила, действующая на 
К, 

площадку поперечного сечения, определяется соотношением [13]: О=2л | о,КаК, где о, — физическая 
0 


компонента тензора напряжений. Поскольку начальная деформация предполагается в виде (1), то тензор 
напряжения Коши Т и мера деформации Фингера Г соосны. Причем тензор напряжений — постоянный: 


2 г: ы . г 
Т = а" (а (ее, +ее,)+4°ъ)-с(а“(е,е, +е,е,) +4“) +сС'(е,е, +ее, + ьь,)) : 


Торцы пластины свободны от нагрузки, поэтому о, =0. Отсюда следует условие, связывающее 


коэффициенты закона Фингера 
а 
с® с +-——=0. (10) 
а? 


Принимая во внимание закон состояния, можно записать конкретное выражение условия (10), которое 
устанавливает связь между аи 4. Так, для материала Мурнагана получаем условие: 


ре +[еи- Зи —4у, + зи - 








3 15 
3).—2 + “. + —- +, + (2^-Зу, —4у, _ + (у +у, — а“ =0. 
В частности, если у, =у, =у, =0 ‚то имеем 
4? Ри. 
А+2и 
При рассмотрении закона состояния Блейтца и Ко 4 выражается формулой 
ей 


Таким образом, задавая значения модулей упругости, нули функции Бесселя, из (9) находим 
бифуркационные значения параметра начальной деформации а,. 


Отметим, величина Л =|-а, для сравнительно тонких дисков — малый параметр. Так, в классической 
теории выпучивания пластин критическое значение величины а, определяется формулой 
2 
в (3,8317) [ й 
12+) (и, 


где у — коэффициент Пуассона, а 3,8317 — первый корень функции Бесселя .Л(х). 





* 


При решении начальной краевой задачи предполагается возможность применения метода возмущений 
Синьорини. При этом считаем, что коэффициенты оператора ® зависят от малого параметра Л степенным 
образом. Это означает, что условия на границе для компонент добавочного перемещения задаются Л матрицей 
А= А(Л). В данной задаче имеет физический смысл частичная проблема собственных значений — 
определение низших собственных чисел. Старшие степени Л, которые содержат элементы матрицы А, 
незначительно сказываются на значениях наименьших корней векового уравнения (9). Если при решении 
начальной задачи ограничиться линейной теорией, то А(Л) — регулярный двучлен [16]. 

Численные эксперименты при исследовании устойчивости нелинейно-упругих тел не слишком 
большой относительной толщины подтверждают этот вывод. Следовательно, характеристическое уравнение с 
нелинейным вхождением параметра (9) можно заменить характеристическим уравнением линейного оператора. 
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Можно применить итерационные процессы, сходящиеся к одному собственному числу, где в качестве нулевого 
приближения выбиралось число, близкое к значению в теории пластин. 

Обсуждение и заключение. С помощью полученных выше уравнений выполнен численный анализ 
влияния физической и геометрической нелинейности на величину верхнего критического значения параметра 
а,. Расчеты реализованы в среде МаЙаЬ для различных вариантов задания модулей упругости, относительной 


толщины пластины и числа волнообразования п. Обнаружено, что во всех случаях потери устойчивости диска 


минимальному критическому значению параметра а отвечает первая форма выпучивания, которая 


соответствует корню функции Бесселя 3,8317 [17]. 
В табл. 1 приведены критические значения параметра а, для пластин относительной толщины 0,05; 0,1; 
0,15; 0,2; 0,25 и 0,3. 




















Таблица 1 
Критические значения параметра а, для пластин различной относительной толщины 
в у Относительная толщина пластины 
и 0,05 0,1 0,15 0,2 0,25 0,3 
| 0,9922 0,9869 0,9800 0,9660 0,9530 0,9330 
0,9972 0,9867 0,9735 0,9549 0,9300 0,8914 
2 0,9977 0,9908 0,9801 0,9664 0,9504 0,9332 
0,9967 0,9869 0,9690 0,9410 — —- 
3 0,9976 0,9907 0,9800 0,9662 0,9502 0,9330 
0,9965 0,9856 0,9665 0,9372 — — 
я 0,9977 0,9907 0,9798 0,9659 0,9497 0,9324 
0,9881 — — 0,9965 0,9613 0,9861 
5 0,9914 0,9865 0,9783 0,9668 0,9520 0,9350 
0,9929 0,9729 0,9341 0,9114 —- — 
6 0,9978 0,9912 0,9800 0,9675 0,9520 0,9340 
0,9985 0,9941 0,9871 0,9775 0,9651 0,9539 
7 0,9985 0,9941 0,9859 0,9735 0,9560 0,9320 
0,9985 0,9941 0,9863 0,9749 0,9594 0,9396 





























Здесь цифрами обозначены перечисленные ниже материалы. 

Сталь Вех 535 (\.=1,09, и=0,818, у =-1,75 , у, =-2,40 ‚ у, =-1 69). 
Сталь 50ХГСМ2Ф (^. = 1,129, и = 0,803, у, =-2,8 , у, =-2,7, у, =-187). 
Сталь Неба 37 (^. =1,11, и=0,821, у, =-3,58, у, =-2,82, у,=-177). 
Сталь Неа АТУ (^=0,87, и=0,716, у, =0,34, у, =-5,52, у, =-Ь0). 
Бронза бериллиевая (^,= 1,042, и=0,49, у, =-4,0, у, =-17, у, =-0,6). 
Оргстекло (^=0,39, и=0,186, у, =-0,078, ъ, =-0,07, ъ, =0,047). 


О о 


В первых шести вариантах рассмотрен закон состояния Мурнагана. Последний вариант представляет 
результаты для материала Блейтца и Ко. Верхнее число в ячейке таблицы относится к случаю, в котором не 
учитывается физическая нелинейность, т.е. у =у, =у, =0 в законе Мурнагана; р =0 в законе Блейтца и Ко. 


Второе (нижнее) число учитывает физическую нелинейность. В законе Мурнагана были использованы данные 


дин 


> 12 
для модулей упругости третьего порядка из [14] в единицах 10 о 


В последней строке табл. 1 приведены результаты для закона Блейтца и Ко, который выбран в 
упрощенном варианте (уравнение Ноулса — Стернберга): коэффициент Пуассона принимался равным 0,25, а 
уточняющий модуль В = 0,45. Прочерк означает, что критическое значение а, не обнаружено. 

Анализ результатов позволяет сделать некоторые выводы. При малых относительных толщинах диска 
точная теория и линейная теория пластин дают одинаковые критические значения параметра а.. 
Геометрическая нелинейность заметно влияет при относительных толщинах больше 0,1. Еще сильнее 
сказывается физическая нелинейность. Однако следует с осторожностью выбирать значения упругих модулей 
третьего порядка. Например, в четвертом варианте (сталь Нес!а АТУ) и для тонких пластин не были 
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обнаружены надежные критические значения параметра а,, хотя для тех же значений модулей упругости Ляме 


наблюдаются приемлемые значения, свидетельствующие о потере устойчивости. Как и в задаче устойчивости 
нелинейно-упругой сферы из материала Блейтца и Ко [15], константа В слабо влияет на величину критического 


значения а,. 
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